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Abstract. Dimuon production has been studied in a series of fixed-target experiments at Fermilab 
during the last two decades. Highlights from these experiments, together with recent results from 
the Fermilab E866 experiment, are presented. Future prospects for studying the parton distributions 
in the nucleons and nuclei using dimuon production are also discussed. 
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INTRODUCTION 

During the last two decades, a series of fixed-target dimuon production experiments 
(E772, E789, E866) have been carried out using 800 GeV/c proton beam at Fermilab. 
At 800 GeV/c, the dimuon data contain Drell- Yan continuum as well as quarkonium 
productions (SP¥, *F', and T resonances). The Drell- Yan process and quarkonium pro- 
ductions often provide complementary information, since Drell- Yan is an electromag- 
netic process involving quark-antiquark annihilation while the quarkonium production 
is a strong interaction process dominated by gluon-gluon fusion at this beam energy. 

The Fermilab dimuon experiments covers a broad range of physics topics. The Drell- 
Yan data have provided information on the antiquark distributions in the nucleons |Q1 
2] and nuclei O, 3]. These results showed the surprising results that the antiquark 
distributions in the nuclei are not enhanced contrary to the expectation of models 

which predict nuclear enhancement of meson clouds. Later, the measurement of the 
Drell- Yan cross section ratios p + dj p + p clearly established the flavor asymmetry 
of the d and u distributions in the proton, and the ^-dependence of this asymmetry 
was determined Q2Q. Pronounced nuclear dependences of quarkonium productions were 
observed for J/^F, and T resonances p, |6|]. The nuclear dependence of Drell- Yan 
cross sections has also provided information on the energy loss of quarks traversin g th e 



nucleus 04LI7J]. In addition, the decay angular distributions for Drell- Yan BEI9D, J/ 1 !' ttlOfl, 
and T resonances ifllll have been measured. 

In this article we first discuss the subject of the flavor structure of the sea quark distri- 
butions in the nucleons. The observation of a striking flavor asymmetry of the nucleon 



sea, inspired by the seminal work of Tony Thomas II 1211 . has profound implications on 
our current knowledge on the parton substructures in the nucleons. Some recent results 
from the dimuon production experiments are then presented. Finally, prospects for future 
experiments to study flavor structures of the nucleons and nuclei will be discussed. 



FLAVOR STRUCTURE OF LIGHT-QUARK SEA 



The earliest parton models assumed that the proton sea was flavor symmetric, even 
though the valence quark distributions are clearly flavor asymmetric. The similar masses 
for the up and down quarks suggest that the nucleon sea should be nearly up-down 
symmetric. The issue of the equality of u and d was first encountered in the Gottfried 
integral fl 1 311 , given as 

h = J [F/(x)-F 2 M (x)] J xdx =\ + \j Q [u p (x)-d p (x)]dx, (1) 

where F£ and F£ are the proton and neutron structure functions. Under the assumption 
of a u, d flavor-symmetric sea in the nucleon, the Gottfried Sum Rule [13], Iq = 1/3, is 
obtained. The most accurate measurement 111411 of the Gottfried integral gives 0.235 ± 
0.026, significantly below 1/3. This surprising result has generated much interest. It 
is remarkable that, already in 1983, Tony Thomas predicted a large excess of d to u 
as a direct evidence for a pionic component in the nucleon. In his paper uM, Tony 
commented on "the lack of experimental information about the shapes of s(x), s(x) and 
d(x) — m(jc)", and concluded his paper "with a plea for better measurements of these 
three quantities in the free proton". 

The shape of the d{x) — u(x) was later measured in proton-induced Drell-Yan and 
semi-inclusive DIS experiments. At Fermilab, the E866/NuSea Q2Q Collaboration mea- 
sured the DY cross section ratios for p + d to that of p + p: 

o DY (p + d)/2o DY (p + p) ~ [l+d(x)/u(x)] jl. (2) 

This ratio was found to be significantly different from unity for 0.015 < x < 0.35, 
showing an excess of d to u over an appreciable range in x. 

Many theoretical models, including meson-cloud model, chiral-quark model, Pauli- 
blocking model, instanton model, chiral-quark soliton model, and statistical model, have 
been proposed to explain the d/u asymmetry. Details of these various models can be 



found in some review articles 1151 116Q. These models also have specific predictions for 



the spin structure of the nucleon sea 111711 . In the meson-cloud model, for example, a 
quark would undergo a spin flip upon an emission of a pseudoscalar meson 
n°(uu,dd) + w^, — > n + (ud) +d^, — > K + + s^, etc.). The antiquarks {u,d,s) are 
unpolarized (Au = Ad = As — 0) since they reside in spin-0 mesons. The strange quarks 
(s), on the other hand, would have a negative polarization since the up valence quarks in 
the proton are positively polarized and the — » K + process would lead to an excess 
of s^. By considering a vector meson (p) cloud, non-zero u,d sea quark polarizations 



with Ad — Au > were predicted 118 



The Pauli-blocking model Q20Q predicts that an excess of q^(q^) valence quarks would 
inhibit the creation of a pair of q^q^ (q^q^) sea quarks. Since the polarization of the 
u(d) valence quarks is positive (negative), this model predicts a positive (negative) 
polarization for the u(d) sea {Au > > Ad). 



In the instanton model [21], the quark sea can result from a scattering of a valence 
quark off a nonperturbative vacuum fluctuation of the gluon field, instanton. The cor- 
relation between the sea quark helicity and the valence quark helicity in the 't Hooft 



effective lagrangian (i.e. ifl leads to a d^) naturally predicts a positively (negatively) po- 
larized u(d) sea. In particular, this model predicts Il22n a large Aw, Ad flavor asymmetry 
with Aw > Ad, namely, / [Aii(x) — Ad(x)]dx — | f [d(x) — w(jc)] 

In the chiral-quark soliton model [2J, llU], the polarized isovector distributions 
Aii(x) — Ad(x) appears in leading-order (N 2 ) in a l/N c expansion, while the un- 
polarized isovector distributions u(x) — d(x) appear in next-to-leading order (N c ). 
Therefore, this model predicts a large flavor asymmetry for the polarized sea 
[Au(x) — Ad(x)] > [d{x] -u(x)}. 

In the statistical model [12511 . the momentum distributions for quarks and antiquarks 
follow a Fermi-Dirac distributions function characterized by a common temperature and 
a chemical potential /! which depends on the flavor and helicity of the quarks. It can be 
shown that /i^ = — and fl^ = —fJLqf. Together with the constraints of the valence 
quark sum rules and inputs from polarized DIS experiments, this model leads to the 
prediction that d > u and Au > > Ad. 

Measurements of Au{x) and Ad(x) are clearly of great current interest. Both the 
HERMES S and the COMPASS |27|] collaborations have reported results on the 
extraction of Au(x) and Ad(x) from polarized semi-inclusive DIS data. These results 
show that Aw, Ad are small, but with large uncertainties. A recent global analysis ^2$\ 
of polarized DIS and polarized p — p interaction indicates that Au(x) > 0, Ad(x) < 0, 
and | Aw(jc) | < | Ad(x) \ . This interesting result suggests that the sea-quark polarization is 
flavor-asymmetric and of opposite sign compared to the unpolarized case. Additional 
data are anticipated for W -boson production at RHIC [29]. The longitudinal single-spin 
asymmetry for W production in polarized p + p collision is sensitive to Au(x) and Ad(x) . 

While various theoretical models can describe the general trend of the dj u asymmetry, 
they all have difficulties [30, HI] explaining the Fermilab E866 data at large x (x > 0.2), 
where dj u drops below 1 . However, the E866 large-* data suffer from large statistical 
uncertainties, and more precise measurements are needed. The 120 GeV Main Injector 
at Fermilab and the new 30-50 GeV proton accelerator, J-PARC, present opportunities 
for extending the d/u measurement to larger x (0.25 < x < 0.7). For given values of x\ 
and X2 the DY cross section is proportional to l/s, hence the DY cross section at these 
lower energies are significantly larger than at 800 GeV. A definitive measurement of 
the d/u over the region 0.25 < x < 0.7 could be obtained for an upcoming experiment 
E906 Il32ll at Fermilab and a proposed measurement [3311 at J-PARC. 

To disentangle the d / u asymmetry from the possible charge-symmetry violation ef- 
fect [|34f| . one could consider W boson production in p + p collision at RHIC. The ratio 
of the p + p — > W + +X and p + p — >■ W~ +X cross sections is sensitive to d/ u. An im- 
portant feature of the W production asymmetry in p + p collision is that it is completely 
free from the assumption of charge symmetry [35fl . Another advantage is that it is free 
from any nuclear effects. Moreover, the W production is sensitive to dj u flavor asymme- 
try at a Q 2 scale of ~ 6500 GeV 2 /c 2 , significantly larger than all existing measurements. 
This offers the opportunity to examine the QCD evolution of the sea-quark flavor asym- 
metry. A recent study showed that W asymmetry measurements at RHIC could provide 
an independent determination of d/u pal . 

While it is generally assumed that the gluon distributions in the proton and neutron 
are identical, this assumption has not been tested experimentally. A possible mechanism 



for generating different gluon distributions in the proton and neutron, as pointed out by 
Piller and Thomas [37|], is the violation of charge symmetry in the parton distributions in 
the nucleons [35]. Unlike the electromagnetic Drell-Yan process, quarkonium produc- 
tion is a strong interaction dominated by the subprocess of gluon-gluon fusion at 800 
GeV beam energy. Therefore, the T production ratio, o(p + d — >T)/o(p + p — > Y), is 
expected to probe the gluon content in the neutron relative to that in the proton [|38|]. 
The o(p + d)/2(j (p + p) ratios for T production with 800 GeV proton beam have been 
reported recently [381], and they are consistent with unity, in striking contrast to the 
corresponding values for the Drell-Yan process. The T data indicate that the gluon dis- 
tributions in the proton and neutron are very similar. These results are consistent with no 
charge symmetry breaking effect in the gluon distributions. 



TRANSVERSE SPIN AND DRELL-YAN PROCESS 

The study of the transverse momentum dependent (TMD) parton distributions of the 
nucleon has received much attention in recent years as it provides new perspectives on 
the hadron structure and QCD B9I1 . These novel TMDs can be extracted from semi- 
inclusive deep-inelastic scattering (SIDIS) experiments. Recent measurements of the 



SIDIS by the HERMES [40] and COMPASS [41|] collaborations have shown clear 



evidence for the existence of the T-odd Sivers functions. These data also allow the first 



determination 114211 of the magnitude and flavor structure of the Sivers functions and the 
nucleon trans versity distributions. 

The TMD and transversity parton distributions can also be probed in Drell-Yan ex- 
periments. As pointed out [|43j] long time ago, the double transverse spin asymmetry 
in polarized Drell-Yan, Ajt, is proportional to the product of transversity distributions, 
hi(x q )hi(xq). The single transverse spin asymmetry, A^, is sensitive to the Sivers func- 
tion Il44ll . fij{x) of the polarized proton (beam or target). Even unpolarized Drell-Yan 
experiments can be used to probe the TMD distribution function, since the cos20 az- 
imuthal angular dependence is proportional to the product of two Boer-Mulders func- 



tions [45], h\ {x\)h\ [x2). A unique feature of the Drell-Yan process is that, unlike the 
SIDIS, no fragmentation functions are involved. Therefore, the Drell-Yan process pro- 
vides an entirely independent technique for measuring the TMD functions. Furthermore, 
the proton-induced Drell-Yan process is sensitive to the sea-quark TMDs and can lead to 
flavor separation of TMDs when combined with the SIDIS data. Finally, the intriguing 
prediction [46|1 that the T-odd TMDs extracted from DIS will have a sign-change for the 
Drell-Yan process remains to be tested experimentally. 

No polarized Drell-Yan experiments have yet been performed. However, some infor- 
mation on the Boer-Mulders functions have been extracted recently from the azimuthal 
angular distributions in the unpolarized Drell-Yan process. The general expression for 
the Drell-Yan angular distribution is Il47n 

da „ ? v ? 

— oc l+Acos 2 + usin20cos<i> + -sin z ecos2<i>, (3) 
dil 2 



where and are the polar and azimuthal decay angle of the l + in the dilepton rest 
frame. Boer showed that the cos 20 term is proportional to the convolution of the quark 



and antiquark Boer-Mulders functions in the projectile and target [48]. This can be 
understood by noting that the Drell-Yan cross section depends on the transverse spins 
of the annihilating quark and antiquark. Therefore, a correlation between the transverse 
spin and the transverse momentum of the quark, as represented by the Boer-Mulders 
function, would lead to a preferred transverse momentum direction. 

Pronounced cos 20 dependences were indeed observed in the NA10 rf49ll and 
E615 Il50n pion-induced Drell-Yan experiments, and attributed to the Boer-Mulders 
function. The first measurement of the cos 20 dependence of the proton-induced Drell- 
Yan process was recently reported for p + p and p + d interactions at 800 GeV/c |@]. 
In contrast to pion-induced Drell-Yan, significantly smaller (but non-zero) cos20 az- 
imuthal angular dependence was observed in the p + p and p + d reactions. While 
the pion-induced Drell-Yan process is dominated by annihilation between a valence 
antiquark in the pion and a valence quark in the nucleon, the proton-induced Drell-Yan 
process involves a valence quark in the proton annihilating with a sea antiquark in the 
nucleon. Therefore, the p + p and p + d results suggest [@, SI] that the Boer-Mulders 
functions for sea antiquarks are significantly smaller than those for valence quarks. 



FUTURE PROSPECTS AT FERMILAB AND J-PARC 

Future fixed-target dimuon experiments have been proposed at the 120 GeV Fermilab 
Main Injector and the 50 GeV J-PARC facilities. As discussed earlier, the Fermilab E906 
experiment will extend the dj u asymmetry measurement to larger x region. Another goal 
of this experiment is to determine the antiquark distributions in nuclei at large x using 
nuclear targets. New information on the quark energy loss in nuclei is also expected. An 
advantage of lower beam energies is that a much more sensitive study of the partonic 
energy loss in nuclei could be carried out using the Drell-Yan nuclear dependence [0]. 

With the possibility to accelerate polarized proton beams at J-PARC [1520 , the spin 
structure of the proton can also be investigated with the proposed dimuon experiments. 
In particular, polarized Drell-Yan process with polarized beam and/or polarized target at 
J-PARC would allow a unique program on spin physics complementary to polarized 
DIS experiments and the RHIC-Spin programs. Specific physics topics include the 
measurements of T-odd Boer-Mulders distribution function in unpolarized Drell-Yan, 
the extraction of T-odd Sivers distribution functions in singly transversely polarized 
Drell-Yan, the helicity distribution of antiqaurks in doubly longitudinally polarized 
Drell-Yan, and the transversity distribution in doubly transversely polarized Drell-Yan. 
It is worth noting that polarized Drell-Yan is one of the major physics program at 
the GSI Polarized Antiproton Experiment (PAX). The RHIC-Spin program will likely 
provide the first results on polarized Drell-Yan. However, the high luminosity and the 
broad kinematic coverage for the large-x region at J-PARC would allow some unique 
measurements to be performed in the J-PARC dimuon experiments. 
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